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We theoretically study the quantized anomalous Hall effect (QAHE) in skyrmion crystal (SkX) without ex¬ 
ternal magnetic field. The emergent magnetic field in SkX could be gigantic as much as ~ 4000T when its 
lattice constant is ~ Inm. The band structure is not flat but has a finite gap in the low electron-density regime. 
We also study the conditions to realize the QAHE for the skyrmion size, carrier density, disorder strength and 
temperature. Comparing the SkX and the system under the corresponding uniform magnetic field, the former 
is more fragile against the temperature compared with the latter since the gap is reduced by a factor of ~ 1/5, 
while they are almost equally robust against the disorder. Therefore, it is expected that the QAHE of SkX system 
is realized even with strong disorder at room temperature when the electron density of the order of one per a 
skyrmion. 


Magnetic skyrmion is a topological spin texture in 
ferromagnet^. After the early theoretical proposals in mag- 
netJ^, the study of magnetic skyrmion is growing rapidly 
since it was discovered experimentalljH^. The periodic ar¬ 
ray of skyrmions, i.e., a skyrmion crystal (SkX), is real¬ 
ized at interfaceJ^ or in bulk chiral magnets such as B20 
compoundJ^. An emergent magnetic field, generated in the 
background skyrmion spin texture. Namely, a skyrmion has 
one magnetic flux $o = h/e acting on the conduction elec¬ 
trons coupled to it. When the skyrmions form a periodic lat¬ 
tice, i.e., a SkX, the emergent magnetic field reaches ~ 4000T 
assuming the uniform averaged flux for the skyrmion size and 
the lattice constant of SkX of the order of ~ Inm. The ef¬ 
fective magnetic field is proportional to A“^, where A is the 
skyrmion radius. Since the size of the skyrmion is 3nm for 
MnGd33 18nm for MnSP, and 70nm for FeG^^, the corre¬ 
sponding emergent magnetic field is ~ HOOT, 28T, and IT, 
respectively. 

This emergent magnetic field leads to the Hall effecP. Most 
of the studies focus on the Hall effect in metallic systems with 
large electron densit}®!!. This so-called topological Hall 
conductivity axy is usually small compared with the longi¬ 
tudinal conductivity i-C-. the Hall angle (Jxyl<Jxx is typi¬ 
cally of the order of 10“^ at mosPl 

Up to now, we regard skyrmions as the source of the real 
space emergent magnetic field. When the size of the skyrmion 
becomes comparable to the mean free path, it is expected that 
the crossover from the real to momentum space Berry curva¬ 
ture occurs. One can regard the case of pyrochlore ferromag- 
net as the limit of large emergent magnetic field, where the 
spin chirality is defined for each unit cell of tetrahedrorP^. In 
this case, there is no Landau Level (LL) formation, but the 
band structure is formed by taking into account the solid an¬ 
gle of the spin, and the intrinsic anomalous Hall effect appears 
whose conductance is given by the integral of the Berry cur¬ 
vature in momentum spac^i^. A more drastic example is the 
quantized anomalous Hall effect (QAHL) in magnetic topo¬ 
logical insulator (TI), where the surface state with gap open¬ 
ing due to the exchange coupling to the magnetic ions pro¬ 
duces the quantized Hall conductance of without the ex¬ 
ternal magnetic fielcP. 

It is expected that the SkX offers an ideal laboratory to 
study QAHL from a unified viewpoint since one can change 


the size of the skyrmion, the mean free path, and even the 
carrier concentration by gating at the interface, to reveal the 
crossover between real and momentum Berry curvature and 
stability of the quantized Hall conductance as these conditions 
are changed. 

In this paper we theoretically explore the emergence of the 
QAHL in the SkX without external magnetic field. Landau 
Levels are not formed in the present system since the emer¬ 
gent magnetic field is not uniform. Nevertheless, the band 
structure contains several well separated bands in the low 
electron-density regime, where each band has a Chern num¬ 
ber C = — 1. Consequently the emergence of the QAHL is 
predicted. We point out that the lowest and next-lowest bands 
are well described by the Dirac theory. On the other hand. 
Hall plateaux disappear in the large electron-density regime 
because of the overlap of bands. We also clarify the condi¬ 
tions of QAHL for the skyrmion size, carrier density, disorder 
strength and temperature. 

The model: We start with a free electron system 
coupled with t he background spin texture Ui by Hund’s 
couplin^i^E3ill The Hamiltonian is given by the double¬ 
exchange model, 

H = — J^^riiclcrci, (1) 

ij 

where c| (ci) is the two-component (spin up and spin down) 
creation (annihilation) operator at the i site, is the trans¬ 
fer integral between nearest-neighbor sites, J is the Hund’s 
coupling strength between the electron spin and background 
spin texture and cr denotes the Pauli matrix. When Hund’s 
coupling is strong enough J , the spin of th e hop ping 

electron is forced to align parallel to the spin texturd^^I®. The 
wave-function |x('r)) of the conduction electron at r corre¬ 
sponding to the localized spin n{r) is given by 

|x(r)) = (cos^^,e*‘^(’’)sin^^)*, (2) 

where we have introduced the polar coordinate of the spin 
configuration n = (cos(/)sin0,sin0sin0,cos0). Then the 
effective transfer integral is obtained, 

^eff = ^(XilXi) = cos (3) 








2 


(a) band structure 



(b) band structure 



-4 --- 

r X M r 


(c) band structure (uniform) 
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FIG. 1: The band structure in the presence of SkX. We have set A = 
2. The horizontal axes are momentum and ky, while the vertical 
axis is the energy. There are 4A^ = 16 bands. We show (a) a bird’s 
eye view and (b) a cross section of the lowest 8 bands, (c) The band 
structure of a tight-binding model with the uniform magnetic flux, 
which forms the almost flat Landau levels except near the center of 
the energy 0. 
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where = {(j)i — (j)j){\ — cos ‘ ^ ^ )/2 is the vector potential 
generated by the spin between the i and j sites, and 0^ is 
the angle between the two spins. The effective tight-binding 
Hamiltonian is obtained as 

^3 


where d\ (di) is the spinless creation (annihilation) operator 
at the i site. One can easily see that the spin chirality is absent 
for collinear and coplanar spin alignment. 

Skyrmion crystal: We consider a background spin texture 
n(r) made of a square SkX. Each skyrmion has a nontrivial 
topological numbeiHSI. The skyrmion profile is well assumed 
as 9{r) = 7r(l — r/A) for r < A and 6{r) = 0 for r > A. The 
emergent magnetic fiel d is p roduced by the spin texture since 
it has a finite solid angleH^SEl 


bz{r) 


n 

Ye 


n ■ {Oxfi X dyTi) 


h TT 

-- 7 sin7r(l 

2erX ^ 


-r/A) (5) 


for r < A and = 0 for r > A. It does not depend on the 
azimuthal angle cj). The total magnetic flux is 

d^rb^{r) = $o, (6) 

with = h/e, and is independent of the skyrmion radius A. 

Band structure: The band structure can be obtained by 
numerically diagonalizing the Hamiltonian (|^ in the unit cell 
with the size 2A x 2A in the presence of the SkX. The Bril- 
louin zone is given by —7r/2A < kx < tt/2X and —7r/2A < 
ky < 7r/2A. Here we take the unit with the lattice constant 
a = 1. Then the unit cell has dimensionless area (2A)^. We 
show the band structure in the presence of SkX in Figs|T])a) 
and (b), which are warped and the number of bands is 4A^. 
There are finite gaps between two successive bands for lower 
bands. However, the band overlap starts at a higher band. In 




FIG. 2: (a) Band gaps and overlaps for various sizes A of the 
skyrmions. (b) The (2A)“^ dependence of the lowest band gap A. 
The gap is proportional to (2A)“^. The band gap of the system with 
the corresponding uniform mean magnetic flux is shown in the blue 
line. 


Fi^TJc) shown the band structure of the corresponding uni¬ 
form mean magnetic field. Here, since one skyrmion exists 
per area 4A^, the mean magnetic field b is given by $o/4A^. 
It is seen that there is almost no energy dispersion, i.e., the 
Landau level formation, with the energy separation Aq which 
slightly depends on the Landau level index. 

We show the gaps and overlaps of bands in Figj^a) for 
various skyrmion radius A, where the bands are marked in 
color bar and the band gaps are denoted by the white blanks. 
We find that the band gap shows a scaling behavior for the 
skyrmion radius. The band gap between the lowest and 
second-lowest bands always opens. 

We show the A dependence of the lowest band gap A in 
Fig|^b). It is proportional to 1/A^. There is some devia¬ 
tion from the linear fit in small A, which is probably due to 
the finite size effect. We can compare this band gap with the 
Landau level separation Aq shown in Fig. 1(c). The lowest 
Landau level gap Aq is given by 0.71t for A = 2, while the 
lowest gap of the QAHE in SkX is given by A = O.llt. The 
linear relation in the small 1 /4A^ region of Fig. |^b) indicates 
the relation A = Ao/5. 

Berry curvature: We focus on the lowest and second- 
lowest bands [Fig|^a)]. We show the spin direction deter¬ 
mined by ('0(fc)|tT|'0(fc)) for conduction and valence bands in 
Fig|^c) and (d), respectively, where the anti-vortex structure 
is evident at the M point. It implies that the spin texture has a 
nontrivial Berry curvature. 

We may define a "gauge potential" in the momentum space, 
ak{k) = —i {ipik)] dk ^{k)), for Bloch state |'0(fc)), which 
is properly called the Berry connection. Then we may de¬ 
fine the "magnetic field" or the Berry curvature by bz{k) = 
dk^o,y{k) — dk^ax{k). The Chem number is the integral 
of the Berry curvature over the first Brillouin zone, C = 
^Jd‘^kb,{k). 

We show the momentum dependence of the Berry curvature 
in Fig|^ It takes a large value in the vicinity of the M point, 
which is {kx,ky) = (7r/2A, 7r/2A). The sign of the Berry 
curvature is negative for the lowest band, while it is positive 
for the second-lowest band. Accordingly, the sign of the Berry 
curvatures are opposite between the lowest and second-lowest 
bands. 

Each band has one unit Chem number C = — 1. This fact 
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(a) band structure 


-2 



(b) Berry curvature (c) spin direction 



FIG. 3: (a) Band structure of the lowest and second-lowest bands 
in the vicinity of the M point. It has the Dirac-cone shape, (b) 
The momentum distribution of the Berry curvature in the lowest and 
second-lowest bands. The Berry curvature takes the largest value at 
the M point, where the Ghent numbers C = —1/2 and -1-1/2 are 
generated for these bands, respectively. The residual Ghent numbers 
C = —1/2 and —3/2 arise away from the M point. The spin di¬ 
rections of the band structure of (c) the conduction and (d) valence 
bands. They clearly exhibit the anti-vortex structure at the M point. 


can be understood in terms of adiabatic connection from the 
Landau level. In our system, the emergent magnetic field has 
a space dependence. We consider an adiabatic pass from the 
uniform magnetic field to the space-dependent magnetic field 
induced by SkX. This deformation has no singularity. Thus 
the band structure in Fig. 1 can be obtained from the defor¬ 
mation of the Landau level. Each Landau level has the same 
Chern number. Accordingly, each band of our system also has 
the same Chern number, which is C = — 1. 

The structure of the lowest and second-lowest bands has 
precisely the shape of the Dirac cone in the vicinity of the 
M point {kx,ky) = (7r/2A, 7r/2A) as in Fig[^a). This sug¬ 
gests that electrons are described by the Dirac theory H = 
hv{—kxO'x + kyay) + maz, where the mass m is related to 
the gap A between the lowest and second-lowest bands by 
A = 2 |to|. 

Hall cond uctan ce: The conductance is calculated by the 
Kubo formul^SEll 
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where n and m are the band indices and f{x) is the Fermi 
distribution function. We show the Hall conductance in Fig|^ 
calculated using the Kubo formula (j^. At the zero tempera¬ 
ture, it is quantized to be the Chern number when the Fermi 

energy is inside the gap, a^y = \ Lnifiiied where is 
the Chern number of the nth band. Below the band gap, the 
Hall conductance decreases, while it increases above the band 
gap. This is due to the fact that the sign of the Berry curvature 
is opposite between the two adjacent bands. The total Chern 
number can be very large for large A, which is distinct from 



FIG. 4: (a) The Hall conductance as a function of the chemical poten¬ 
tial. The horizontal axis is the chemical potential, while the vertical 
axis is the Hall conductance a^y The Hall conductance is quantized 
as marked by arrows. There is a peculiar structure of dips in the con¬ 
ductance at the values where it is quantized, (b) Zoom up of (a) near 
the band edge. It is well fitted by the dotted curve obtained in the 
Dirac theory. 



FIG. 5: (a) The Hall conductance at finite temperature. Though it is 
not quantized, the peculiar dip structure is observable even at room 
temperature, (b) Zoom up of (a) near the band edge. 


the QAHE in magnetic topological insulator^® 

This behavior of the Hall conductance across the gap can 
be interpreted by the Dirac theory, which gives 


_ / —1/2 for IpI < \m\ 

-m/{2\n\) for jp] > \m\ ' 


( 8 ) 


It describes the peculiar behavior of the Hall conductance 
quite well as in Fi^b), where the conductance is quantized 
inside the gap and continuously changes outside the gap show¬ 
ing a dip structure. 

Finite temperature: We show the Hall conductance at fi¬ 
nite temperature in Fig|^ It is evident that the temperature 
scale is given by the gap A in Fig. 1(b), which is around O.lf 
for A = 2. Although the quantization of plateau is broken at 
ksT = O.Olf ~ O.OleL ~ lOOK, the peculiar dip struc¬ 
ture toward the formation of plateau is clearly visible even at 
ksT ~ AOOK. Consequently, our prediction of QHE with¬ 
out external magnetic field can be experimentally observable 
at room temperature. 

Disorder effect: Up to now, we focus on the pure system, 
and the Berry curvature in momentum space. Now we intro¬ 
duce the disorder potential as given by 


Ffijnp ^ ^ Ui dj, 


I 


(9) 
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FIG. 6: (a) The Hall conductance as a function of the chemical poten¬ 
tial for various disorder strength. The horizontal axis is the energy, 
while the vertical axis is the Hall conductance a^y Disorder strength 
is indicated by color. Hall plateaux disappear as disorder increases. 
However, the Hall plateau at a^y = —e^jhis transformed to the 
Hall plateau at a^y = 0. (b) The Hall conductance in the vicinity of 
the band edge. 


where Ui takes a uniform random distribution for—V < Ui < 
V. We have calculated the Hall conductivity axy for the sys¬ 
tem of size 8x8, including 4 unit cells. Figure 6 shows the 

numerical results, which clearly shows the plateau transition 

2 

occurs at V =t for the lowest band, i.e., (j^y = ~ x • ^hat 

when |(t| > 5 • IT’ the two-parameter scaling trajectory con¬ 
verges to the quantized Hall state. Taking this criterion, the 
lowest occupied band, which is most promising candidate to 
observe the QHE in SkX, can support the plateau up to V =t. 
This magnitude of the random potential is similar to the sep¬ 
aration between LLs in Fig. 1(c) and is much larger than the 
gap A in Fig. 1(b) which is = 0.It in the present case. This is 
because the dispersion gives the stability or robustness of the 
extended state. For higher energy bands, the Hall conductiv¬ 
ity is reduced more slowly as V increases, and the criterion 
\a\ > 5 • X satisfied up to 1/ = 3t as shown in Fig.6(a). We 
also calculated quantity Pn for the eigenfunction (ilV'n) with 
energy £„ defined by 

Pn=Y.\^\^r.)\\ ( 10 ) 

i 

which measures the extent of the wavefunction. Namely, 


when {i\4’n) extends over the M sites, Pn ^ M There- 
fore, the localization length is estimated by P„ in 2D. We 
show in the inset of Fig. 6 the averaged ' for the lowest 
band (purple curve) and all the bands (black curve) as a func¬ 
tion of V. It seen that the disappearance of the QHE roughly 
corresponds to the localization of the wavefunctions less than 
4, where 4 is the size of the unit cell (skyrmion). Therefore, 
since the mean free path £ is always longer than the local¬ 
ization length, the QAHE is observed only in the case where 
£ > 2A, i.e., the Berry curvature in momentum space is more 
appropriate picture rather than that in the real space. 

Discussions: We have shown that the QAHE occurs in 
metallic SkX due to the emergent magnetic field induced by 
the spin texture. It is shown that there is a gap between the 
lowest and second-lowest bands. The electron density should 
be of the order of one electron per one skyrmion, whose elec¬ 
tron density is given by 1 /4A^. In order to realize experimen¬ 
tal situations, it is necessary to reduce the number of electrons. 
In this sense, an interface of dilute magnetic semiconductors 
will be promising. The emergent magnetic field and the asso¬ 
ciated energy gap is proportional to 1/A^. Namely, skyrmions 
with smaller radius are better for realizing QAHE. 

The QAHE is robust against disorder as in the case of sys¬ 
tem with corresponding uniform magnetic field with Landau 
levels. Considering the very large mean magnetic field of 
~ 4000T for the skrmion size A ~ Inm, the disorder does not 
destroy the QAHE so seriously. This is because the change 
in the Chern numbers occurs only via the pair annihilation, 
and the appreciable energy dispersion of the band in SkX pro¬ 
tects the Chem number. On the other hand, the band gap A of 
SkX is smaller than Aq, i.e., that of the corresponding system 
with uniform magnetic field, and hence the stability against 
the thermal excitation is reduced. A is estimated as Ao/5, but 
as mentioned above, the emergent magnetic flux induced by 
SkX is gigantic and even the QAHE even at room tempera¬ 
ture is expected. 
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